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Abstract 

Extending the gauge-invariance principle for r functions of the standard bi- 
linear formalism to the supersymmetric case, we define M = 1 supersymmet- 
ric Hirota operators. Using them, we bilinearize SUSY nonlinear evolution 
equations. The super-soliton solutions and extension to SUSY sine-Gordon 
are also discussed. As a quite strange paradox it is shown that the Lax in- 
tegrable SUSY KdV of Manin-Radul-Mathieu equation does not possess N 
super-soliton solution for N > 3 for arbitrary parameters. Only for a partic- 
ular choice of them the N super-soliton solution exists. 

I. INTRODUCTION 

Supersymmetric integrable systems constitute a very interesting subject and as a conse- 
quence a number of well known integrable equations have been generalized into supersym- 
metric (SUSY) context. We mention the SUSY versions of sine-Gordon p|, [EJ, KP-hierarchy 
||, KdV Q], J4|, Boussinesq |j] etc. We also point out that there are many generalizations 
related to the number M of fermionic independent variables. In this paper we are dealing 
with the Af = 1 superspace. 
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So far, many of the tools used in standard theory have been extended to this frame- 
work, such as Backlund transformations ||, prolongation theory, hamiltonian formalism 
grasmmannian description ||], r functions ||, Darboux transformations jlQfl . The phys- 



ical interest in the study of these systems have been launched by the seminal paper of 
Alvarez-Gaume et. al [|TT|] about the partition function and super- Virasoro constraints of 2D 



quantum supergravity. Although the r function theory in the context of SUSY pseudodiffer- 
ential operators was given for the SUSY KP-hierarchy |§, the bilinear formalism for SUSY 
equations was very little investigated. We mention here the algebraic approach using the 
representation theory of affine Lie super- algebras in the papers of Kac and van der Leur |12| . 



Kac and Medina [13[ the super-conformal field theoretic approach of LeClair |14| . Anyway 



in these articles the bilinear hierarchies are not related to the SUSY hierarchies of nonlinear 
equations. 



This paper which is an extended version of we consider a direct approach to SUSY 
equations in a M = 1 superspace rather than hierarchies namely extending the gauge- 
invariance principle of r functions for classical Hirota operators. Our result generalize the 



Grammaticos-Ramani-Hietarinta |T7[ theorem, to SUSY case and we find TV = 1 SUSY Hi- 
rota bilinear operators. With these operators one can obtain SUSY-bilinear forms for SUSY 
KdV equation of Mathieu Q and also it allows bilinear forms for certain SUSY extensions of 
Sawada-Kotera-Ramani Hirota-Satsuma [19], KdV-B ||, Burgers, mKdV and M = 2 



-superspace SUSY sine-Gordon equations p5|] . Also the gauge- invariance principle allows 
to study the SUSY multisoliton solutions as exponentials of linears. A very interesting fact 
which happens is that SUSY KdV equation of Mathieu does not have 3 supersoliton solution 
for arbitrary choice of solitary waves although it possesses Lax pair J|. Only for special 
combination of parameters the equation admits N soliton solutions. Although this seems to 



be a quite strange paradox, Liu and Manas [10], found also super-soliton solutions for SUSY 



KdV equation in terms of pfaffians only for certain wave parameters. This fact shows that 
Hirota integrability and Lax integrability are different in the SUSY context. 

The paper is organized as follows. In section II the standard bilinear formalism is briefly 



discussed. In section III supersymmetric versions for nonlinear evolution equations are 
presented and in section IV we introduce the super-bilinear formalism. In the last section 
we shall present the bilinear form for SUSY KdV-type equations, super-soliton solutions and 
several comments about extension to M = 2 SUSY sine-Gordon equation. 



II. STANDARD BILINEAR FORMALISM 



The Hirota bilinear operators were introduced as an antisymmetric extension of the usual 



derivative [plj] , because of their usefulness for the computation of multisoliton solution of 
nonlinear evolution equations. The bilinear operator Y) x = d xi — d X2 , acts on a pair of 
functions (the so-called "dot product") antisymmetrically: 

Dz/ • g = (d Xl - d X2 )f(x 1 )f(x 2 )\ xl = x ^ x = f'g - fg'. (2.1) 

The Hirota bilinear formalism has been instrumental in the derivation of the multisoliton 
solutions of (integrable) nonlinear evolution equations. The first step in the application is a 
dependent variable transformation which converts the nonlinear equation into a quadratic 
form. This quadratic form turns out to have the same structure as the dispersion relation 
of the linearized nonlinear equation, although there is no deep reason for that. This is best 
understood if we consider an example. Starting from paradigmatic KdV equation 

u t + Quu x + u xxx = 0, (2.2) 

we introduce the substitution u = 2d x log F and obtain after one integration: 

F x tF F x Ff + F XXXX F AlF xxx F x + SF XX = 0, (2-3) 

which can be written in the following condensed form: 

(D.D, + B 4 X )F • F = 0. (2.4) 

The power of the bilinear formalism lies in the fact that for multisoliton solution F's are 
polynomials of exponentials. Moreover it displays also the interaction (phase-shifts) between 
solitons. In the case of KdV equation the multisoliton solution has the following form: 



N 

F = expd^^i + ZjAi^i), (2.5) 

fi=0,l i=l i<j 

where rji = kiX — kft+r]!- and expAij = ( ^'~^ J ) 2 which is the phase-shift from the interaction 
of the soliton "z" with the soliton "j". 

This picture can be generalized to any bilinear equation of the form 

P(D S )F • F = (2.6) 

where P is any polynomial and x = (t, x, y, ...). The 1 soliton solution is F = 1 + e v , where 
7] = kx + const. This solution holds if 

P(k) = 0. 

This is a condition on the parameters k of n and is called the dispersion relation. If the 
parameter space is n-dimensional then the above equation defines an n — 1 dimensional 
submanifold called dispersion manifold. 
Hirota ansatz for 2 soliton solution is 

F = 1 + e m +e m + A 12 e Vl+m 

where r/s are defined before and A 12 are function of k\ and A; 2 each one giving the coordinates 
of some point in the dispersion manifold. Substituting the ansatz in the equation (|2l|) and 
taking into account the dispersion relation we find 

A 12 P{k 1 + k 2 ) + P{k x - k 2 ) = 0. (2.7) 

Generally speaking, the great majority of bilinear equations possess 2 soliton solution but 
only integrable equations possess 3 or more soliton solutions. The form of the N soliton 
solution is, 

N 

F = exp /liTji + AjViVj) 

H=0,l i=l i<j 



More about the bilinear equations are in |22 



4 



A very important observation (which motivated the present paper) is the relation of 
the physical field u = 1d 2 x log F of KdV equation with the Hirota function F: the gauge- 
transformation F — > e P x + wt F leaves u invariant. This is a general property of all bilinear 
equations. Moreover, one can define the Hirota operators using the requirement of gauge- 
invariance. Let's introduce a general bilinear expression, 

N 

A N (f,9) = j:^~ l fW x g) (2.8) 
and ask to be invariant under the gauge-trasformation: 

A N (e 9 f,e 9 g)=e 29 A N {f,g) 9 = kx + out + ...{linear s). (2.9) 



Then we have the following, [I7| 

Theorem: A^(f,g) is gauge-invariant if and only if A^(f, g) = Y) x f»g i.e. 

(n) 

Ci = co(-l) 1 

V 1 ) 

and Co is a constant and the brakets represent binomial coefficient. 

We must point out that in the whole paper the natural number N (which will denote 
number of solitons, number of terms, exponents etc.) is different from M which is related 
to supersymmetry or superspace. 



III. SUPERSYMMETRY 

The supersymmetric extension of a nonlinear evolution equation (KdV for instance) refers 
to a system of coupled equations for a bosonic u(t,x) and a fermionic field £(£, x) which 
reduces to the initial equation in the limit where the fermionic field is zero (bosonic limit). 
In the classical context, a fermionic field is described by an anticommuting function with 
values in an infinitely generated Grassmann algebra. However, supersymmetry is not just a 
coupling of a bosonic field to a fermionic field. It implies a transformation (supersymmetry 
invariance) relating these two fields which leaves the system invariant. In order to have a 
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mathematical formulation of these concepts we have to extend the classical space (x, t) to 
a larger space (superspace) (t, x, 9) where 9 is a Grassmann variable and also to extend the 
pair of fields (u, £) to a larger fermionic or bosonic superfield $>(t,x,9). In order to have 
nontrivial extension for KdV we choose $ to be fermionic, having the expansion 

$(t,x,6)=£{t,x) + 6u(t,x). (3.1) 

The Af = 1 SUSY means that we have only one Grassmann variable 9 and we consider 
only space supersymmetry invariance namely x — > x — X9 and 9 — > 9 + A (A is an anti- 
commuting parameter). This transformation is generated by the operator Q = dg — 9d x , 
which anticommutes with the covariant derivative D = dg + 9d x (Notice also that D 2 = d x ). 
Expressions written in terms of the covariant derivative and the superfield $ are manifestly 
supersymmetric invariant. Using the superspace formalism one can construct different su- 
persymmetric extension of nonlinear equations. Thus the integrable (in the sense of Lax 
pair) variant of TV = 1 SUSY KdV is [| 

$ t + D 6 $ + 3L> 2 ($D$) = 0, (3.2) 

which on the components has the form 

u t = -u xxx - Quu x + 3f f xx 

£t = ~£,xxx ~ 3£a;W ~ 3£,U X . (3-3) 

Another integrable variant of SUSY KdV equation which is very important in applications 
to supersymmetric matrix models is SUSY KdV-B equation ||, namely 

$ t + D 6 <5> + 6D 2 $D$ = 0, (3.4) 

leads to a somewhat trivial system in which the fermionic fields decouple from the bosonic 
equation which reduces then to the usual KdV. 

We shall discuss also the following supersymmetric equations, although we do not know 
if it is completely integrable in the sense of Lax pair ($ is also a fermionic superfield). 
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• J\f = 1 SUSY Sawada-Kotera-Ramani, 

$t + L> 10 $ + D 2 (10D$D 4 ® + 5D 5 $$ + 15(D$) 2 $) = 0. (3.5) 

• M = 1 SUSY Hirota-Satsuma (shallow water wave) 

D A <$> t + $ tJ D 3 $ + 2D 2 $D$ t - L> 2 $ - $ t = (3.6) 

• JV = 1 SUSY Burgers 

$t + + $ ra = (3.7) 

A very important equation from the physical consideration is the SUSY sine-Gordon. 
We are going to consider the version studied by Kulish and Tsyplyaev |2f| . There are other 



integrable versions of SUSY sine-Gordon emerged from algebraic procedures |24[ . In this case 
one needs two Grassmann variables 9 a with a = 1, 2 and the supersymmetry transformation 
is 

a> = - i\^e, e' a = e a + \ a , v = 1,2. 

Here, X a is the anticommuting spinor parameter of the transformation and A = (A 1 , A 2 ), 
A a = \ l 3(i(T2) l3a , 7° = icT2, 7 1 = o"i, 7 5 = 7°7 1 = cr 3 . We use the metric g^ v = diag(—l, 1) 
and (7j are the Pauli matrices. The superfield has the following expansion: 

Q ) = + z^(^) + (3.8) 

where and F are real bosonic (even) scalar fields and ip a is a Majorana spinor field. The 
SUSY sine-Gordon equation is: 

DD$ = 2i sin $, (3.9) 

where -D a = c^a + i('~f' J '6) a d f j i and on the components it has the form: 

(7^ + cos0)V> = O 
(foxx - <j>tt = 2 ( sin ( 2< /0 _ #^ sin 0) • (3-10) 
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IV. SUPER-HIROTA OPERATORS 



In order to apply the bilinear formalism on these equations one has to define a SUSY 
bilinear operator. We are going to consider the following general A/" = 1 SUSY bilinear 
expression 



N 



S N (f,g)=Y.c l (D N ~ l f)(D*g), 

i=0 



(4-1) 



for any N, where D is the covariant derivative and /, g are Grassmann valued functions 
(odd or even). We shall prove the following 

Theorem: The general M = 1 SUSY bilinear expression ( |4.1|) is super-gauge invariant 
i.e. for G = kx + uot + linears (£ is a Grassmann parameter) 



S N (e e f,e e g)=e^S N (f,g) 



20 , 



if and only if 



»(<+!) 



N 



ct = c (-1) j|/I 
where the super-binomial coefficients are defined by: 

if (JV,i) ^ (0,l)mod2 



iV 



' [N/2] ^ 



V [i/2] j 







otherwise 



|/| is the Grassmann parity of the function / defined by: 



l/l 



1 if / is odd (fermionic) 



if / is even (bosonic) 

and [k] is the integer part of the real number k ([k] < k < [k] + 1) 

Proof: First we are going to consider N even and we shall take it on the form iV = 2P. 
In this case we have: 

N P P-l 

S N (f,g) =Y,c.(D N -*f)(D*g) = J2c 2l (d p -*fWg) + £ c 2 j+i(d p ~ j ~ 1 Df)(d j Dg) 

i=l i=0 j=0 
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Imposing the super-gauge invariance and expanding the covariant derivatives we obtain: 



EE 



E C 2* 



P—n—m 



k^ 1 



(d m f){d n g) + 



n>0m>0 ^1=0 




v m ) 


) 




+ E E A 


( P-l 


fj ) 


f p-j~\\ 




\ 


E C 2i+1 


iP—n'—m' 


-i 


n'>0 m'>0 


\r 


\ n 'i 






/ 



+ E E a(-i) 



/i+i 



n/>0 m'>0 



(P-l t \ 

E C 2i+1 
i=0 



n' 



V 



771 





(p-1 






\ 


+ EE 


E 


hi 




ti>Om>0 








) 



p-j-A 



(d m f)(d n Dg)+ 



P—n'—m' — l 



(d m Df)(d n g) + 



J 



( 



\ 



m 



,P—n—m—l 



J 



(d m Df)(d n Dg) 



J 



= J2c 2i (d p - i f)(d i g) + j:c 2j ^(d p -^Df)(&iDg) 

i=0 j=0 

where A = £ + 9k. From this, we must have for every m, n subjected toO<n<i<P — m 
and j < P — w! . 









(p-i\ 


p 






E C 2i 








i=0 






\ m ) 



k 



P—n—m 



C2n$P-n- 



(4.2) 



Also due to the fact that the supergauge invariance has to be obeyed for every / and g we 
must have C2j+i = The discrete equation (|4.2|) was solved in JL7|. Its general solution is 
given by: 



C2i = C (-l) 1 



c 2j+l — 



(4.3) 



In the case oi N = 2P + 1 we proceed in a similar manner and we obtain the following 
system: 



p 


(') 


(p-i\ 


E C 2i 






i=0 


[nj 


\ m J 



k 



P—n—m 



C2n$P-n- 



(4.4) 



j=0 



(A 

[n) 



( P-j-l\ 



m 



P—n—m~l 



—n—m—1 



l) l/l c 2i + c 2i+1 = 



(4.5) 



(4.6) 



This system has the following solution: 



c 2i = co(-iy 



\ l J 



c*+i = c (-l) i+1+ l^l 



\ l J 



(4.7) 



The relations (|4.3| ), (|4.7| ) can be written in a compact form as 



Cj = c (-l 



A 7 



and the theorem is proved. We mention that the super-bilinear operator proposed by 
McArthur and Yung [TjJ is a particular case of the above super-Hirota operator. 
We shall note the bilinear operator as 



S N (f,g) := S^f.g 

In the Appendix 1 we list several simple properties of this super-Hirota operator. 

V. BILINEAR SUSY KDV-TYPE EQUATIONS 

SUSY KdV of Mamn-Radul-Mathieu 

In order to use the super-bilinear operators defined above we shall consider the following 
nonlinear substitution for the superfield: 



&(t,x,0) = 2D 3 \ogr(t,x, 



(5.1) 
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where r is an even superfield. Introducing in SUSY KdV ( |3.2|) and integrating with respect 
to x we obtain the following 

2Dd t log r + 3{ (2D 3 log r) (2d 2 x log r) } + 2D 7 log r = 



Using the properties fl6.4|), fl6.5p, ( |6.6| ) we find 



2 I *y 2 



2^(— ) = o 

r 



(5.2) 



Now we are using the following property of the classical Hirota operator |L9 



3 a = D 3 x a»b _ D g a«6Dg6«6 



In our case a = Dr and b = r. Accordingly 



V r 



Plugging into (|5.2j) we find the following super-bilinear form 



(S s D t + S 7 x )r . r = 0, 



(5.3) 



which is equivalent with the form found by McArthur and Yung W5 



S x (D t + T>l)T*T = 0. 



(5.4) 



In order to find the super-soliton solutions we are going to use the classical perturbative 
method namely the series 



r = l +e/ (D +e 2 / (2) +e 3 / (3) 



(5.5) 



where /w are even functions. Equating the power of e we find: 



for e, 



D(d t + d?)fW = 



(5.6) 
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• for e 2 , 



2D(d t + dl)f^ = -S s (D t + D^)/« . /«. (5.7) 



for e 



3 



■ 



+ <9 3 )/< 3 ) = -S«(D t + D2)/« • / (2) . (5.8) 

• for e 4 , 

2D(d t + dl)f^ = -2S X (D, + D»)/« • / (3) - S X (D, + D 3 )/( 2 ) . /( 2 > (5.9) 
and so on. 

Now if we take = ^-kH+e^) the equation flyj) is satisfied, / (2) = 0, / (3) = 0.... and 



the series (|5.5|) truncates. So, the 1 supersoliton solution is given by 
for every h and £. 

Introducing in the super-bilinear equation the 1 soliton solution F = 1 + 

exp (fcx + cut + (9) one finds the dispersion supermanifold equation: 

P(fc,u;,C) = (C + ^)( W + fc 3 ) = 

which imposes to = — k 3 for every £. 

In order to find 2 super-soliton solution we take = e' ?1 +e' ?2 where ^ = kiX — kft + 8Q 
The equation ( |5.7| ) becomes 

2L>(d t + dl)f^ = Ghhih - fc 2 )[(Ci - fa) + - A^e^" 2 (5.11) 

Taking into account that r is a Grassmann even function, must be also even. Ac- 
cordingly, the general solution of (|5.11|) has the form 



f {2) = [mi2{h, k 2 , Ci, 6) + 0n 12 (h, k 2 , (i, Ca) 

12 



e m+m 



mi2 and h\2 being Grassmann valued even and odd functions. Due to the fact that we have 
two Grassmann parameters (£1 and (2) the most general expressions for m\ 2 and n u are 
given by, 

mi 2 = (kt - hf/ih + k 2 f + /c 2 )Ci6 

n 12 = a (k lt k 2 )Ci + Kh, k 2 )( 2 

The above forms for m 12 and n 12 could also be obtained by expanding in power series of Ci 
and ( 2 and taking into account that in the bosonic limit (Q — > 0), m\ 2 — > (ki — k 2 ) 2 / (ki + k 2 ) 2 
(ordinary KdV interaction term) and h\ 2 — > 
Introducing in the equation (|5.11|) we find 



/ 



(2) 



h - k 2 



o h-fa ; ; na (ki - k 2 ){k 1 ( 2 - k 2 ( 1/ 
(fci + k 2 y 



(h + k 2 ) 



h + k 2 J 

Introducing the above forms for /W and in ( |5.8|) one obtains /^ 3 ^ = and then 
= 0, = 0... i.e. the series truncates. So the 2 super-soliton solution is given by; 



r ( 2 ) = 1 + e Vl + e m + A 12 e 



(5.12) 



where 



A 



12 



'Asi — Aia^ l0 ^i — A: 2 44 „(&i - k 2 )(k 1 ( 2 - k 2 (i) 

+ . , no S1S,2 + — 



(A;i + fc 2 ) 



(A;i + h) 



(5.13) 



One can easily see that the classical general procedure for finding the interaction term given 
by the equation ( |2.7|) 

A 12 P(h + k 2 ) + P{k x -k 2 ) = 

does not work in the SUSY case. 

In order to find 3 supersoliton solution we consider 

/« = e m + e m + e V3 



and the equation for /( 3 ) (|5.8|)becomes: 



D{d t + dl)f^ 



S.Dt + SDie" 1 • A 23 e M + e ^ • A 13 e r?1+7 ' 3 + e" 3 • A 12 e 



(5.14) 
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The solution /( 3 ) of this equation (which is very complicated) does not cancel the right hand 
side of the equation ( |5.9|) . So, is not zero and the series ( |5.5|) cannot be truncated. 
Accordingly the SUSY KdV equation does not have 3 super-soliton solution in the standard 
form for arbitrary k^s and (Vs. This seems a quite strange paradox, because SUSY KdV is 
integrable in the sense of Lax. 

Anyway, imposing the constraint kiQ = kjQ for every i and j it is easy to prove (see 
Appendix 2) that SUSY KdV possesess the following N-soliton solution 

N 

T (N) = exp^/i^ + ^Ay/i^), (5.15) 

//=0,1 i=l i<j 

where 

rji = k t x - k^t + 9Q + rjP 

exp ^- = (trt) 

"'iCj kjCi 



Solutions with constraints on parameters have been found also by Liu and Manas [IU|, using 
SUSY Darboux transformation. 

SUSY KdV-B 

This situation is completely different in the SUSY KdV-B case 

$t + £> 6 $ + 6D 2 $D$ = 0. (5.16) 
With the same nonlinear substitution 

$ = 2D 3 log F 

we obtain the ordinary form 

(D t n x + n i x )F*F = o 

which has N-super-soliton solution ( |2.5| ) the fermionic contribution being only an additive 
phase i.e. rji = kiX — kft + Q9. 
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SUSY Sawada-Kotera-Ramani 



For M = 1 SUSY Sawada-Kotera-Ramani (|3.5|) using the nonlinear substitution, 

$ = 2D 3 log r(t,x,9) 



we shall find the following super-bilinear form: 



;S x D t + Si 1 )r.r = 



(5.17) 



In a similar way we find the 2 super-soliton solution 



-( 2 ) = 1 + e m + e m + 



k\ — k 2 \ 2 k\ — k\k 2 + k 2 ki — k 2 k'f — k\k 2 + k. 



+ 2- 



k 1 + k 2 J kl + k^ + kl (fci + k 2 ) 2 k\ + k x k 2 + k 



rCiC 2 



x(l+2 /f 1 ~f lC2 )e^ 
fci - k 2 

Also with the same constraint kiQ = kjQ we find 



N 



T = J2 ex P <X ^ + Z AijWi)' 

H=0,1 i=l i<j 



where 



rji = k iX - k\t + 6Ci + vf ] 



exp Ai 



k{ — kj \ k^ kikj -\- kj 
ki -\- kj J kf -\- k{kj -\- k 2 



(5.18) 



SUSY Hirota-Satsuma 



For Af = 1 SUSY Hirota-Satsuma equation using the nonlinear substitution: 



$ = 2Dlogr(t,x,9) 



one obtains the super-bilinear form: 



'S 5 r D 4 - Sl - S x D t )r • r = 



(5.19) 



The 2 supersoliton solution is given by: 



r ( 2 ) = 1 + e" 1 + e" 2 + 



ki — k 2 \ . , k\ — k 2 , , ; j 

" > Mi 2 + 2 /; : ; 2 NO M 12 CiC 2 



fci + fc 2 



(fci + ^) 5 



'1 + ofl fc 2Cl fclC2 x „ 1+r;2 
4 A*- A* j 
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where 

(fci - k 2 f + hhjjh - k 2 f - (kf - i){kj - 1)] 

12 (h - k 2 f - kM(ki - k 2 f - (k\ - \)(k\ - 1)] 
and r\i = kix — k^t/ (kf — 1) + Q9. With the constraints k { (j = kjQ the equation admits also 
N soliton solution with 

SUSY Burgers 

An interesting case is the SUSY extension of the Burgers equation, namely 

$ t + $.0$, + = 

It is welknown that the classical Burgers equation can be linearized via Cole-Hopf transform. 
If we try the nonlinear substitution (which is the natural supersymmetrization of the Cole- 
Hopf transform): 

$ = 2D log r(t, x, 9) 



we find: 



S x B t T • t + 2B 2 x Dt • r = 0. 



This form is not super-gauge invariant. Accordingly we are forced to use two functions for 
substitution, namely 

^_ G(t,x,0) 
F(t,x,9) 

where G is an odd Grassmann function and F is an even one. Using the relation (|6.7|) from 
Appendix 1 we obtain the following gauge-invariant super-bilinear form 

(D t + D 2 X )G*F = 

B 2 X F •F = S 3 X G»F 
This system admits the following 1 super-shock solution: 

G = 2(C + 9k)e i - kx - k2t+ t e \ F = 1 + e {kx - kH+ ^ 
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We can ask ourselves if it is possible to obtain super-bilinear forms for SUSY equations 
of the nonlinear Klein-Gordon type. In fact the SUSY sine-Gordon equation( |3.9| ) can be 
written in the following form: 

[D T ,D x MT,X,9,9 t ) = 2ism$(T,X,9,9 t ) (5.20) 

where we have introduced the light-cone variables X := i[t — x)/2, T := i(t + x)/2, and 
9 := 9%, 9 2 := —9 t . Co variant derivatives are Dx '■= dg + 9dx, := dg t + 9 t dr and the 
square braket means the commutator. Using the nonlinear substitution (G and F are even 
functions) 

$ = 2* log (J), 
we find the following quadrilinear expression 

2i{F 2 (G[D T , D X ]G - [D T G, D X G}) - G 2 (F[D T , D X ]F - [D T F, D X F})} = F 4 — G A 

It is easy to see that the bilinear operator 

S X tT mr:=r[D T , D x ]r - [D T r, D x r] 

is super-gauge invariant with respect to the super-gauge 

g 6 ._ e (kx+ivt+8(+e t n+liniars) 

Accordingly we can choose the following super-bilinear form, formally the same with stan- 
dard sine- Gordon equation, 

SxtG*G=^(F 2 -G 2 ) 

SxtF • F = ^(G 2 - F 2 ) (5.21) 

but, it is not clear how to compute the super-kink solutions. 

iFiom these examples it seems that gauge-invariance is a useful concept for bilinear 
formalism in the supersymmetric case, though there is no deep reason for that. As a conse- 
quence we was able to bilinearize two supersymmetric equations of KdV type and the SUSY 
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sine-Gordon. The case of SUSY versions for mKdV, NLS, KP etc. requires further investi- 
gation because it seems that only certain supersymmetric extensions are super-bilinearizable. 
Although we do not know if the SUSY extension of Sawada-Kotera proposed above are inte- 
grable in the sense of Lax, it admits super-bilinear form and only 2 super-soliton solution for 
arbitrary choice of solitary waves. The strange fact is that SUSY KDV equation of Mathieu 
which is known to be Lax integrable does not admit Ngeq3 supersoliton solution in the 
canonical form. Probably a singularity analysis implemented on the super-bilinear form will 
reveal the connection between Hirota-integrability and Lax-integrability. 

VI. APPENDIX 1 

In this section we are going to list several properties of the super-Hirota bilinear operator 
which are useful in deriving bilinear forms. 

STf*9 = B^f.g (6.1) 

S!f +V 1 • e" 2 = [Ci - C2 + 0(h - k 2 )\ (h - k 2 ) V' 1+ " 2 (6.2) 

S™ +ii . e" 1 = (-l) N+1 (( + 6k)k N e^ = (-l) N+1 S 2 x N+1 ei • 1 (6.3) 
where rji = kiX + 9Q and Q are odd Grassmann numbers. 

2Dlogr = — (6.4) 

T 



2 J D 3 logr=Sl^ (6.5) 



2Z^logr= SJ3 ; 2 r,r (6.6) 

where r is an even Grassmann function Moreover if G and F are Grassmann functions (with 
F even) then 

o/Gx S?.G • F \ r \G &..F * F . . 

D 3 ( f ) = ^p,— - {-^p^p,— (6-7) 
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VII. APPENDIX 2 



In this section we are going to sketch the proof of the formula for N supersoliton solution 
for bilinear SUSY KdV equation. We are rely on the proof of Hirota for ordinary N-soliton 
solution in the case of KdV |[21|| , ||23|| . Thus, introducing the expression of the N supersoliton 
solution 

TV 

T (N) = exp(J2^iVi + J2 A v^j)' 

/i=0,l i=l i<j 

r)i = k t x - kft + 9Q + 



into the super bilinear form 



(S x B t + S 7 x )r . r = 0, 



and taking into account the properties ( |6.1D and (|6^ ) we find 

TV TV TV TV 

E E {(£(* - ^)a,)((-) E(^ - + (Eto - ^)Ai)(E(w - ^)*i) 3 }x 

/i=0,l i/=0,l i=l i=l i=l i=l 

TV 

x exp (E()"i + ^t)^ + EO^i + ^My) = 

i=l i<j 

where A = ^ + Since /ij, z/j = 0, 1 it is clear that we only have exponential terms of the 
form 

n m 

exp (^77;+ 2r]i),0<n<N,n<m<N 

i=l i=n+l 

Next we show that the coefficient of this general exponential term is zero, the coefficient 
being given by, |^TJ p3j : 



TV TV TV TV n 

A = E {(-E^)(E^ 3 ) + (E^ A i)(E^) 3 }Il(^-^%) 2 

o=0,l i=l i=l 8=1 8=1 i<j 

where — ^ — z/j. We do not go into details because the procedure of deriving the above 
coefficient goes absolutely in the same way as in the case of ordinary KdV equation. This 
is due to the fact that the interaction term is the same for KdV and SUSY KdV. 
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Because A = Q + 9ki the coefficient A becomes on the components 



N 



N 



N 



N 



a = a + a 1= {(-E^)(E^ 3 ) + (E^)(E^^) 3 }Il(^^-^^) 2 + ( 7 - 2 ) 



(T=0,1 



i=l 



i=l 



i<j 



N 



N 



N 



N 



+0 E {(-E^)(E^ fc3 ) + (E^)(E^) 3 }ri(^^-^^) 2 

(T=0,1 i=l i=l i=l i=l i<j 

Hirota proved that the second term is zero [2~T| , The first term is also zero because, 

using the property Ci&j = Cj^i, it can be written as 



A = = 
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